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Abstract 

We study the local conserved charges in integrable spin chains of the XYZ type with 
nontrivial boundary conditions. The general structure of these charges consists of a bulk 
part, whose density is identical to that of a periodic chain, and a boundary part. In contrast 
with the periodic case, only the charges corresponding to the interactions of even number 
of spins exist for the open chain. Hence, there are half as many charges in the open case as 
in the closed case. For the open spin-1/2 XY chain, we derive the explicit expressions of all 
the charges. For the open spin-1/2 XXX chain, several lowest order charges are presented 
and a general method of obtaining the boundary terms is indicated. In contrast with the 
closed case, the XXX charges cannot be described in terms of a Catalan tree pattern. 
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1. Introduction 

Recently we have obtained tlie structure of tlie local conserved charges for the XXX 
chain [1] and some of its generalizations: the isotropic su{N) spin chain and the octonionic 
model [2] . Explicit expressions - in the form of a Catalan tree pattern - have been derived 
for all the charges. Different extensions of these results can be considered. In the present 
work we study the modification of the Catalan tree pattern for open (finite) chains.0 

As for closed chains, there is a transfer matrix formalism for open chains from which 
the conserved quantities can be obtained by power expansion in terms of the spectral 
parameter. For closed chains, this approach is not an effective way of deriving the explicit 
expressions of the conserved charges. In the open case, the situation is even worse as 
the open transfer matrix is, roughly speaking, the trace of the square of the closed chain 
monodromy matrix (up to boundary terms). A somewhat surprising consequence of this 
squaring process is that half of the local conserved charges of the closed chain disappear 
when the chain is cut open. More precisely, as we will show below, the expression for 
the open conserved charges contain a bulk part and a boundary part. In the infinite chain 
limit, the boundary part is irrelevant: the bulk part is thus necessarily the full closed chain 
charge expression. Hence, a charge exists in the open case only if suitable boundary terms 
can be added to the bulk part to preserve the commutativity of the total charge with the 
hamiltonian. 

It turns out that conservation preserving boundary terms do not always exist. For 
instance, for the XY closed chain, there are two infinite sequences of local charges, that 
is, {Hn^} and {Hn^}, for all integer values of n, where n indicates the highest number 
of adjacent spin factors of the charge density. In the open case, there are no appropriate 
boundary counterterms for half of these charges. There remains only one charge for each 
value of n. For the closed XXX chain, we know that there exists one local charge for 
each positive integer n. However, for the open XXX case, only the charges with n even 
survive. In particular, does not exist! 

A major technical difficulty encountered in deriving explicit expressions for the con- 
served charges for open chains is caused by the non-existence of a boost operator, whose 
commutation with the hamiltonian would generate the higher order charges. As a re- 
sult, there is no systematic and effective way of obtaining the expressions for the first few 

^ As far as the pattern of charges is concerned, the infinite open chain is not distinguishable 
from the closed case. 
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charges. In spite of this, one could try to find them by brute force and see whether an 
underlying pattern emerges. From the expressions of H4 and Hq obtained for the XXX 
open chain, this would seem to be so: their boundary parts can be written in terms of 
the same spin polynomials as their bulk part and they can actually be expressed in terms 
of lower order 'bulk charges' localized at the boundaries. However, this is not a generic 
feature: Hg cannot be written in this polynomial basis. 

The article is organized as follows. In the next section, we review the transfer matrix 
formalism. This is used to show the absence of odd degree local conserved charges. A 
second argument is worked out in the following section, from a computational approach: 
the charge is written as a sum of a (known) bulk piece and a set of boundary counterterms 
to be fixed. We show that in the counterterms must cancel a Catalan-tree bulk-type 
charge localized at each boundaries. We then demonstrate that no such terms can be 
found for and then for all H2n+i- The explicit form of the first few nontrivial charges 
is presented. Unfortunately the simple polynomial basis used to describe the closed XXX 
charge is inadequate in the open case, thus making it very hard to find a general expression 
for all the charges. Straightforward generalizations and conclusions are reported in the final 
section. 

2. The general structure of local conservation laws in open chains 

In this section we briefiy recall the transfer matrix formalism and use it to determine 
the general structure of conservation laws in open spin chains. This will be done first for 
the relatively simple spin-1/2 XXX model, and later the extension of the argument to more 
general integrable chains will be indicated. 

The integrable spin chain with nontrivial boundary conditions [3] is described by the 
transfer matrix 



u is the spectral parameter and T{u) is the monodromy matrix, whose trace gives the 
usual (closed) chain transfer matrix: 



(the subscripts 'o' and 'c' refer to open and closed cases respectively). The monodromy 
matrix is constructed from the basic i?-matrix of the model 



to{u) = tr{K+{u)T{u)K_{u)T-\-u)). 



(2.1) 



tc{u) = trT(w) 



(2.2) 



T{u) — Rno • ■ • Rio- 



(2.3) 
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The index i = 1, - ■ • ,N labels a vector space at site number i and refers to internal or 
auxiliary space, over which the trace in (2.1) is taken. K+, K_ are respectively left and 
right boundary matrices. For the XYZ model, the most general form of these matrices 
compatible with integr ability has been derived in [4]. For free open boundaries, to which 
we confine ourselves for the most part of this work, K± = /, the 2x2 identity matrix. 

The transfer matrix can be used as a generating function for the conserved quantities. 
It follows from the Yang-Baxter equation that any derivative of the transfer matrix com- 
mutes with the hamiltonian (which is, by construction, the first logarithmic derivative of 
the transfer matrix, taken at the particular value m = of the spectral parameter). The 
quantities obtained in this way are however nonlocal, i.e., they contain spin interactions 
at distances growing arbitrarily when the length of the chain is increased. To obtain local 
conserved quantities out of the transfer matrix, one needs to take logarithmic derivatives 
[5] . The argument of Luscher, showing the local nature of the logarithmic derivatives of 
tc{u) at tt = 0, remains valid for open chains. However, there are two new features in the 
the open case as compared to the periodic one. First, translational invariance is obviously 
lost. Second, and more surprisingly, all charges with an odd number of interacting spins 
disappear! 

We will first concentrate on the spin- 1/2 XXX chain, whose R matrix takes the simple 

form 

R = ul + P (2.4) 

where P is the permutation operator. The derivatives of the transfer matrix of the closed 
chain are: 



u=0 

C=N- 



where the sum goes over all ordered clusters C of {N — n) points, that is, sets {ii, ■ ■ .ijv-n} 
with ii < 12 < ■ ■ ■ < iN-n- ^fl(C) (and Vl{C) appearing below) denote the cyclic per- 
mutation of the spins of the sequence C to the right (or to the left, respectively). In 
particular, 

tc(0) = 7'fl(A), (2.6) 

where A = {1, 2, . . . , A/"}. The logarithmic derivatives of the transfer matrix yield local 
charges. In particular, using 

tZ\Q) = Vm, (2.7) 



it can be easily seen that 

N N 



^c-'(0)t?)(0) = J2'^Ri{i,i + l}) = J]Pm+i- (2.8) 
The permutation operator above can be written as 

P,,i+1 = |(l + (7i-(7i+i) (2.9) 

(where cr^ denotes a Pauh matrix at site i) and thus (2.8), up to a constant, coincides with 
the hamiltonian of the periodic XXX chain. 

Let denote the (n — l)-th logarithmic derivative of tc{u) at tt = 0: 

<5n=^^lnt,(.x)U^o. (2.10) 

In terms of Vk defined as 

Vk = T-\0)T^''\0) (2.11) 

we have 

= tlVn-l + trpn-livi, . . .,Vn-2), (2.12) 

where Pn-i is a homogeneous polynomial (of degree n — 1) m vi, . . . , Vn-2- These charges 
are translationally invariant and can be written in the form 

N 

Q'n = Y.^n{j) (2.13) 

where the density qn{j) contains interactions of at most n spins on sites {j, j + 1, . . . , j + 
n — 1} (addition understood modulo N). 

Consider now the chain with open free boundaries {K± — I). Conserved charges may 
be obtained from the expansion of the logarithm of the transfer matrix to{u): 

Qn = ^^lnt„(«)U=o. (2.14) 

We will now argue that for n even, the leading term - that with the greatest number of 
interacting spins - in the above expression is the same as in the closed case. Taking the 
trace over the two-dimensional auxiliary space, we obtain 

to(0) = 21 (2.15) 
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where / is the 2^ x 2-^ identity matrix. Therefore, ^(0) is simply 1/2, and 

Qn = ¥ir~'Ho)+Pn-iCA'\o)^ • • • ' ¥^r'Ho)) (2.16) 

where Pn-i is the same polynomial as in (2.12). Consider first n to be even. From (2.1) 
we get 

(0) = 2 tr Vn-i + tr r^-i (2.17) 

where r^-i is a homogeneous polynomial of degree n — 1 in the variables fi, . . . ,fn-2- 
Note that (2.16) and (2.12) both contain the same local term 

tr^,_i ~^P^,({z,z + l,...,z + n-2}) (2.18) 

i 

(with the addition understood modulo for the closed chain). For the closed XXX chain, 
the set of local conserved charges obtained from the transfer matrix has been shown to 
be complete (see [6]). This means that any local conserved quantity in the closed chain 
must be a linear combination of the logarithmic derivatives of the transfer matrix tc{u). 
In particular, this implies that given the leading term in (the one describing n-spin 
interactions), there is a unique set of subleading terms ensuring the commutativity of 
with the hamiltonian. The introduction of nontrivial boundary conditions can be viewed 
as a perturbation of the original periodic chain II and this can only reduce the number 
of conservation laws. The key point is that, as can be seen from (2.16) and (2.12) for n 
even, the leading term of a charge must be the same as for the closed chain. Consider 
now the infinite limit of the open chain. Obviously, in this limit (2.10) must coincide with 
some function of conserved quantities for the infinite periodic chain. But the only local 
conserved quantities in this case are linear combinations of the charges Q^. Therefore, the 
bulk density of the open chain n-spin charge Q° must be the same as the density of Q^. 
It follows that for n = 2m, has the form 

Q2m — Q^2rn + 92m + '?2m5 (2-19) 

where 

N-2m+l 

Qlt^= E ^2m(i), (2.20) 

i=i 

^ If He denotes the hamiltonian of a periodic spin chain with nearest-neighbor interactions, 
He = Xl^i^i'J+i' adding the single perturbation term —hi^N one obtains the open-ended 
chain Ho = He — /ii.w. 
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(i.e., we take only clusters that do not 'jump' over the boundary) and q'2m i^2Tn) stands 
for boundary terms at the left (right) boundary. As we will show in sect. 3, the boundary 
terms involve at most the 2m — 2 spins adjoining the boundary. 

In contrast, when n is odd, i.e., n = 2m + 1, there is no totally local term in Q^: 

Q2m+1 = P2m{'A'^ . . . . . ^.t^--'^) (2.21) 

and the n-spin contribution vanishes. In consequence, the bulk part of Q2m+i zero, 
implying that for n odd, Q° must be trivial. 

The absence of odd-order charges is in fact a simple consequence of the chain reversal 
symmetry of the open chain, i.e., its invariance under the transformation 

O'n CTN-n+l- (2.22) 

This symmetry is specific to the open case: the transfer matrix of the periodic XXX chain 
is not invariant under (2.22): 

tT{RN,o{u)RN-i,o{'^) . . . i?i,o(w)) ^ tv{Ri,o{u)R2,o{u) . . . RnA^))- (2-23) 

Indeed, it is easy to see from the explicit form of the local conserved quantities of the XXX 
chain (see [1]), that the logarithmic derivatives of tc{u) of odd degrees (corresponding to 
charges with an even number of spins) do not change under (2.22), while those of even 
degree (corresponding to odd number of spins) change sign. 

In contrast, the transfer matrix of the open XXX chain is invariant under the operation 
(2.22). Using 

to{u) can be written in the form 

to{u) =tr(i?Ar,o • • • Ri,o{u)R^}i{-u) . . . i^^^-tt) 

(-1)^ (2.25) 

= -^^-^P^^^tr(i?7V,o(w) • • • i?2,o(w)i??,o(w)^2,o(w) • • • i?7V,o(w))- 

Under (2.22), it transforms as 



(-1) 
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to{u) , ir{RiAu) . . . Rn-i,o{u)Rj^ q{u)Rn-i,o{u) . . . Ri,o{u)) (2.26) 

[i -\- u ) 



which, due to the cychc property of the trace, is equal to (2.25). This symmetry excludes 
thus the possibility of building-up the bulk density of Q2m+i from the densities of odd- 
spin charges Q^m^x since their sign is changed under (2.22). Hence, there are no odd-spin 
charges for the open XXX chain. 

The above argument for the spin-1/2 XXX chain may be similarly applied to the open 
XYZ chain. As before, the key step is to observe that, for n even, (2.16) and (2.12) both 
contain the same local leading term with n spins interacting. In the XXX case, one could 
use the completeness property of the system of charges generated by the transfer matrix 
of the periodic chain to show that the bulk density in the open case must coincide with 
the density of the periodic chain. To our knowledge, the completeness of the system of 
charges (2.10) for the general XYZ case has not been established. However, considering the 
infinite limit of the open chain, one may conclude that the bulk density of the conserved 
charge (2.10) must coincide with the density of the corresponding closed chain expression, 
modulo the densities of additional conserved charges (which do not have to be generated 
by the transfer matrix (2.2)). Therefore, the set of conserved quantities generated from 
the transfer matrix still has the general structure (2.19), but the bulk density may now, a 
priori, contain additional contributions. For n odd, there is no n-spin local term in (2.16), 
and in consequence the bulk density of Q2m+i must vanish. As for the XXX case, this can 
be viewed as a consequence of the chain reversal symmetry (2.22). Recall that the basic 
R matrix of the XYZ model satisfies the unitarity requirement 

R{u)R{-u) = p{u)I (2.27) 

where p{u) is some scalar function. Using (2.27), the transfer matrix of the open XYZ 
chain can be rewritten in a form which is manifestly invariant under (2.22): 

to{u) =tr(i?jv,o(«) • . .-R2,o(w)-Ri,o(«)-Rl^o(-^)-^2",o(-'") • ■■^N^oi-'^)) 

1 2 ' ' (2-28) 

= tr(i?i,o(tt) . . . Rn-i,q{u)Rj^ q{u)Rn-i,q{u) . . . i?i,o(«))- 

p[u) 

Note that the above argument does not exclude a possible existence of additional odd-spin 
conserved charges not given by the formula (2.10). Such a situation indeed exists for the 
XY model (see sect. 4). However, this is not expected in the non-degenerate XYZ case, for 
which the system of charges generated from the transfer matrix is likely to be complete. 



3. The calculation of the conserved charges for the open XXX model 

The model-defining hamiltonian for the XXX chain with free open boundaries isil 



N 



H2 = '^(Ji ■ (Ji+i (3.1) 



i=l 

where a stands for the vector (a^, cr^, cr^), cr's being the usual Pauli matrices. As we have 
seen in the previous section, higher-order charges will have a bulk part and a boundary 
part. It is convenient to redefine the charges so that in the bulk they will be described by 
the Catalan tree pattern [1]. (This corresponds to a change of basis from {Qn} to {Hn} 
where such that does not contain any linear combination of lower order Hm^s. This 
change of basis corresponds to taking linear combinations of the logarithmic derivatives of 
the transfer matrix). We thus look for charges in the form 

if. = if,*^"!^ + C + (3.2) 

where and are boundary terms located respectively at the left and right extremity 
of the chains (their precise extension will be evaluated later) ; these are the unknowns to be 
determined. With the above ansatz, we assume the chain to be sufficiently long to prevent 
a mixing of boundary terms from opposite boundaries. Recall that the local integrals of 
motion for the closed XXX model can be expressed as linear combinations of the F^f^^s 
defined as 

K,k= E fr.{C) (3.3) 

Qin,k) g^aj^^fjg for the set of all n ordered sites with k holes (i.e., the sites are not necessarily 
adjacent) and C is some cluster among this set. /n(u, ^2, ■ " " 5 ^n) is defined as 

/o = /i = 0, 

/s = ((^ii X (^ia) ■ (^i3, (3-4) 
fn = {--- {{(^i, X (J,J X ■ ■ •) ■ ai^. 



From now on we omit the superscript 'o' from conserved charges. 
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(Notice that the parentheses are nested toward the left). The closed chain conservation 



laws -ff^'s are 



[n/2]-l k 

k=i e=i 

where the coefficients ak,e are the generalized Catalan numbers 

f 2k-£ \ f 2k-£ \ 

The bulk part of an open-chain conserved charge (of even order) is then given simply 
by restricting the corresponding closed chain expression to those clusters which do not 
'jump over the border'. More precisely, let O^"''^) denote the subset of obtained 
by removing the clusters which overlap the borders. (For instance, for N = 6, = 
(1, 2, 6), (1, 3, 6), (1, 4, 6), (1, 5, 6).) Denoting 

K,k= E /"(^)' (3.7) 

ffMk giygn by formula (3.5) with ^ replaced by j^. 

For the purpose of determining the boundary terms, we can focus on a single extremity, 
or equivalently, consider a semi-infinite chain (N ^ cx), so that i = 1, 2, ■ • ■). The term 
can then be characterized as follows. It is designed to cancel, in the commutator [H2, Hn], 
those terms that would be cancelled for an infinite chain {i G Z) by the contribution of 
the 'link' ctq ■ cti- In the infinite chain commutator [ao ■ 0"i, H^'^], these are precisely those 
resulting terms that live in the semi-infinite chain, i.e., on sites i > 1. Our first step is 
thus to evaluate this commutator: 

[H2, H^^-'^] = -[ao■a^, H^T^] ^ (3.8) 

where TZ^ stands for left remainder. This expression is readily evaluated: it is a sort of 
anti-boost, that is, the inverse of the boost action that generates Hn+i from In the 
boost operation, [7^2, Hn] gives, in addition to Hn+i, some lower order local charges. This 
turns out not to be the case for the anti-boost calculation. We find the remarkably simple 
result 

Wo■a^,Hjr%y,=^iH:L-j^ (3.9) 

L-bdry 

where Hj°^l^ restricted to the left boundary contains only those terms of H^_^l^ whose 
leftmost spin is at site 1. This is proved below. 



9 



Let us first introduce a schematic way of describing the related calculations. To 
represent the monomial fniH,i2, ■ ■ ■ ,in) we use a sequence of dots, with black dots corre- 
sponding to occupied sites. For our purposes, we will only need sequences of dots based 
at site (that is, the first dot represents ii = 0). For instance /6(0, 5) and /4(0, 1, 3, 5) 



are represented respectively by 

/6(0, 5) = (((((70 X CTi) X (72) X 0-3) X 0-4) • 0-5 ~ (3.10) 

and 

/4(0, 1, 3, 5) = (((70 X (7i) X (73) • (75 ~ • • o • o • (3.11) 

The commutation relation 

[(7o • (7i, ((70 X (7i) • (72] = 4z((7i • (72 - (7o • (72) (3.12) 

will then be represented by the following diagram: 

• •m — Aiomm — 4i • o m (3.13) 



The contraction in this diagram indicates taking the value of the commutator of the 'link' 
(i.e., the ao • (T\ operator) with the / polynomial represented by a sequence of dots. In the 
following, we will implement the projection onto the semi-infinite chain by setting to zero 
all spin factors at site z < 0, that is, by droping the terms whose first dot is occupied. This 
projection will be indicated by an arrow. For example, in (3.13) it amounts to setting to 
zero the term ■ a2'- 

• ^ 4i o • • (3.14) 
Here are sample calculations of [gq ■ (7i, iifn"*^] |^^^. For n = 4, 

i?!"^ = i^4,o + i^2,i (3.15) 

It is clear that the only terms in H/^ that will contribute to this commutator, after the 
projection, are those whose leftmost spin is at site 0. The contribution of F4 o is 

• ••• Ai o • • • . (3.16) 
The commutator with ^2,1 does not survive the projection: 

•^•=2z«o« ^0 (3.17) 
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The result is exactly the part of F3 o based at the left boundary, that is, 



(3.18) 

L — bdry 



For n = 5, 



we have 



and 



We thus obtain 



H^^ =F5,o + Fs,i (3.19) 

• • • ^ 4z o • • • • (3.20) 

• • o m — i>4zo«o« (3.21) 

= -MH^^""^^ . (3.22) 

L— bdry 

Having worked out these examples, let us return to the proof of the general identity 
. The contribution of a term F^^k to the commutator is clearly 4i Fn-i^k restricted to 
the left boundary: indeed, only monomials of the form /n(0, 1, ■ ■ ■ , n — 1) contribute in the 
commutator and after commutation and projection, each such term is transformed into 
4i fn-i{l, ■ ■ ■ , n — 1), which builds up Fn-i,k restricted to clusters starting at site 1. (The 
number of holes is unaffected in this process.) If n is odd, the structure of the Catalan 

If n is even, the terms F2 k 



tree is not modified and is mapped to 4i H^l^ 



L— bdry 



do not contribute and we simply strip off the last row of the Catalan tree of to get 

the tree of H^\, restricted to the left boundary. 

Up to this point, we have only found what the boundary terms have to cancel 
when commuted with H2- we require 



L a: TT bulk 



(3.23) 

L — bdry 



The existence of a conserved charge is thus reduced to the existence a suitable satisfying 
this equation. 

In this framework, the non-existence of is simply established: there are no 
build up with 1 or 2 spin factors, that could cancel Ai ai ■ 02 upon commutation with 
c^i ■ c"2 + o"2 ■ c"3- The same argument also implies the non-existence of any odd charge, 
simply because h'^Xx contains a three-spin piece and /i2n+i would have to contain exactly 
the same counterterm required to make a charge Hj,. 



11 



Equation (3.23) implies that the boundary part in may contain only interactions 
of up to n — 2 spins adjoining the boundary. For the spin- 1/2 chain, the most general term 
containing multi-spin interactions is a multilinear polynomial in spin variables, which can 
be equivalently represented in the permutation basis as a linear combination of permu- 
tations of the n — 1 spins at the boundary (with {n — 1)1 arbitrary coefficients). These 
coefficients are then fixed by enforcing the condition (3.23). The first two boundary terms 
are: 

/ii = 4PL(l,2), (3.24) 
=^-Vl{1, 2, 3, 4) - Vl{^, 3, 2, 1) + Vl{1, 2, 4, 3) + Vl{2, 1, 3, 4)) 

(3.25) 

+ 8Pi?(l,3)-8Pfl(2,3). 

A particularly useful set of polynomials in spin variables is provided by the /n(C)'s (where 
C is an ordered cluster (zi < • • ■in)), introduced above as the building blocks of the bulk 
part of Hn- They can be expressed in terms of permutations as follows: 

fn{{ii, 12,..., in}) = 2{-i)^-% . . [[Pi^i^_, , Pi„_.i„_3] • • • , Pi.n] (3.26) 

with 

/2(^l,^2) = 2P^,^, -1. (3.27) 

In particular, 

fsih, i2, h) =2i{VL{ii, «2, ^3) - 'Pl(^3, ^2, h)), 
/4(n,*2,^3,M) = -2{VL{ii,i2,i-i,iA) + 'PL{i4:,i3,i2,il) 

- VL{ii,i2,iA,iz) - VLii3,i4,i2,ii)), 
f6{ii,i2, is, iA, ib,i6) =2(7'l(h, 22, iz, ■is, ■ie) - 'Pl(h, i2, is, ■ie, ■is) 

-^L(n,«2,«3,«5,«6,^4) +^L(n,^2,«3,«6,«5,^4) 

- 'PL{ii,i2, i4, ib, ie, is) + 7^l(h, ^2, ^4, ie, ib, is) (3.28) 
+ VL{ii,i2,i5,i6,i4,is) - ^l(^i, «2, ^5, ^4, ^3) 

- VL{ii,is, i4, id, ie, i2) + VL{ii,is, i^, ie, id, i2) 
-\-'PL{ii,is,ib,iG,iA,i2) - ^L(n, ^3, «6, ^5, ^4, ^2) 
-\-VL{ii,iA,i?„i6,is,i2) - VL{ii,iA,i6,id,is,i2) 

- VL{ii,i5, i&, iA, is, ^2) + VL{ii,i6, ib,iA, is, ^2))- 
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The boundary terms and hg have a simple form when expressed in the fn{C) basis. 
Define 

fn,kij)= E /-(^) (3-29) 

cec("''=)(j) 

where C^'^'^\j) stands for the set of (n, k) clusters based at site j (recall that k indicates 
the number of holes). Then we have: 

hi = 2/2,0(1), (3.30) 

= 2/4,0(1) + 4/2,1(1) - 4/2,0(2). (3.31) 

One might hope that all the higher boundary terms could be similarly expressed in terms 
of the polynomials fn,kij), which would significantly reduce the number of coefficients in 
the ansatz for the boundary part of Hn- It is easily seen to hold for the leading boundary 
term: 

h2m = 2/2m-2,o(l) + lower orders. (3.32) 

Unfortunately, this is not true in general for the subleading terms as can be seen from the 
exact expression for hg 

= 2/6,0(1) - 4/4,0(2) + 2/4,1(1) + 8/2,0(3) + 2/2,0(2) 

(3 33) 

+ 8/2,1(1) + 2/2,2(1) - 8/2,1(2) + 2/4(1, 3, 4, 5) - 2/4(1, 3, 2, 4). 

Since the last term corresponds to a disordered cluster, cannot be rewritten as a linear 
combination of the simple polynomials /n(C) (with C ordered). Higher-order boundary 
terms are rather complicated and no simple pattern seems to emerge. 



4. The XY chain and related models 

In this section, we present the results for the conservation laws of the open chain of 
the XY type. The basic open XY model is defined by the hamiltonian: 

7V-1 

HxY = E + ^y^j^hi- (4-1) 

i=i 

As is well known, by a the Jordan- Wigner transformation, the XY chain (regardless of 
the boundary conditions) can be reduced to a free fermionic theory. Not surprisingly, the 
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conserved charges for free-fermion chains have a particularly simple form. For periodic 
boundary conditions, the conservation laws have been described in [7, 8, 9, 10] (see also 

[I])- 

An interesting feature of the closed XY case is the existence of two independent families 
of conservation laws, which both persist when the model is perturbed by a perpendicular 
(i.e., in the z-direction) magnetic field. As we will see shortly, when the closed chain is cut 
open, only half of each family survives. 

We consider first the particularly simple special case = Ay = 1 (the XX model). 
We introduce the notation (for n > 2) 

efU) = . . . (4.2) 

and 



(4.3) 



(4.4) 



et U) = U) + U) n even, 

= CO')-er(j) ^odd, 

en U) = <^ U) - er U) n even, 

= erO') + er(j) ^odd, 
with 

er(j) = -^|- (4.5) 
For the periodic XX chain, the two mutually commuting families of conserved charges are: 

N 

"Ht = Y.^t{3). (4.6) 
In particular, there are two 'two-spin hamiltonians': 

^^2- = E^K+i-^K+i- (4-8) 

is invariant under global parity transformations, i.e., a" ^ while '^H~ behaves 

like a pseudoscalari. 



'^i^2 is a special case of the Dzyaloshinski-Moriya interaction [11]. 
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The open chain versions of these hamiltonians 

N-l 

(4-9) 

i=i 

no longer commute with each other: 

[H+,H^]^4tat,-4ial (4.10) 

To each hamiltonian, there corresponds a different infinite family of conserved charges. 
Their bulk parts are given by 

N-n+l 

jj±Mik^ E e±(i). (4.11) 

The family of conservation laws for the scalar hamiltonian is: 

= + ht'"- + ht'"^ n even, 

= i/-;^'-^^ + h--^^ + h--^^ n odd, 



(4.12) 



where 



[n/2]-l 

hf^'^^- E ^n-2k{k). (4.13) 

[n/2]-l 

^n^'' = - E <-2k{N-k + l). (4.14) 



For example. 



hp"^ = - (e+(l) + e+(2) + e+(3)) = -(a^alafalcTsVe- + aMalala^Vl (4.15) 
+ a^alala^ + alalalal + a^al + afa|). 

Obviously, /i^'"^ can be obtained from h^'^ via the chain reversal (2.22). The conservation 
of the family H^^ 

[iy+,iyW] = o (4.16) 

can be verified by a straightforward calculation, using 

[Ht, e±(j)] = ±2z (-1)" (ej+,(j - 1) - e^^.ij) + ej_i(j + 1) - ej_i(j)) (4.17) 
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(4.18) 



where 1 < j < N. 

Similarly, one obtains the family of conservation laws for the open chain pseudoscalar 
hamiltonian 

^(2) =H--Mlk + g-;L ^ g-;R ^ ^^^^^ 

^^+;6n//c^^+;L^^+;i^ ^ ^^^^ 

where 

[n/2]-l 

9^''"^= E (-l)'+'et2fcW (4.19) 

k=i 

and g^'^ is obtained from g^'^ by chain reversal. For example, 



^8"' =66(1) + 61(2) + e^iS) = alalalalalal - ala^alalala^ 



(4.20) 



The two families {Hn^} and {Hn^} do not commute with each other. In particular, 
the charges {Hn^} do not commute with H2 ■ Thus only half of the scalar and pseudoscalar 
families of the periodic chain can be modified by adding the boundary terms so that they 
still commute with the scalar hamiltonian of the open chain. Both families are invariant 
under a global spin rotation around the 2;-axis: it is easily checked that all Hn^ 's commute 
with the generator of such rotation, the 2;-component of the total spin, 

TV 

S^ = J2a]. (4.21) 
Therefore, both families survive when the model is perturbed by a magnetic field term 
Consider now the anisotropic case (XY model): 

N-l 

For the periodic XY chain, there are also two distinct families of conservation laws, all 
commuting together. One sequence is given by 

N 

'Hn = J2h^^{j) (4.23) 

with 

hrU) = Ke^U) + \eyy{3) + A,ef_2(j) + A.e^ ^(j). (4.24) 
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This expression gives the bulk density of the conserved charge for the open XY chain 

Hn = if^''^ + h^ + hi (4.25) 

N-n+l 

HT'= E ^n'^ (•?■)• (4-26) 

The border part is given by 



[n/2]-l 

^n=- E {'^^'""'=VMen-2.(fc) + er-2.(/c-l)] 
/c=l 

+ A/'=A,^-^MC-2fe w + e2fc(^ - 1)]} 



(4.27) 



(4.30) 



where 

Pk = 1(1 - (-1)') (4.28) 

with the convention 

C^(0) = 0. (4.29) 

For example 

h^ = - Xyeni) - X.el\l) - A,er(l) - A^ef (1) - A.e^ (2) 
- A,ef (2) - A,er (3) - A.e^ (2) - A.ef (3) - A.ef (2). 

The commutativity of (4.25) with the hamiltonian (4.22) can be established by a direct 
calculation. The second family of conservation laws (containing the pseudoscalar XX 
hamiltonian) can be modified in a similar way to account for the anisotropic deformation. 

In the presence of both the anisotropy and the magnetic field i.e., for the open XYh 
model 

N-l N 

iy = E A.^,VJ+i + Xyay,^, + E ha], (4.31) 

except for the special cases Aa; = Xy (XXh model) and Xy = (the Xh model, equivalent to 
the Ising chain) there exist no counterterms that could be added to the bulk part of H3 to 
account for the open boundary conditions. The first nontrivial charge in the nondegenerate 
case is H4, with 

N-3 N-2 

Hlulk ^ ^ {X,eT{j) + Xyef{j)) + E («3er(i) + /^sef (j)) 

(4.32) 



N-l N 

n 

3 



+ E("2erO-)+/32ef(j)) + «E^. 
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and 

= -XyGlal - Kg\gI - h{\ + + \ylK)(jl, (4.33) 

where 

a^ = -KK^2\y)l\, (4.34) 

a-2 = -2A^ - Xx'^/Xy + K^/K + h'^/Xy + Kt^/h (4.35) 

and /3j = ai{\x ^ Xy) {i = 2, 3), k being arbitrary. 

We end this section with a remark on the open Hubbard model: 

AT-l AT 

H=-2 J2 E[(4.%+i,^ + 4+i,.«i.«)+4t^EKT-V2)(%u-l/2), (4.36) 

where aj ^ and fermionic creation and annihilation operator of an electron of spin 

s at site j, satisfying the anti-commutation relation: 



a 



0'k,s' + Ofc,s'o],s = ^j,fc^s,s', (4.37) 



^3,s = 4 .aj>. (4.38) 



C/ is a coupling constant and 

The integrability of this system has been established by means of the Bethe Ansatz in [12]. 
As is well known, this model may be equivalently regarded as two coupled XX chains: 

N-l N 

(where a and r denote two independent sets of Pauli matrices). The local conserved charges 
of the Hubbard model have been investigated in [1] for the case of periodic boundary 
conditions. For free open boundaries, the first nontrivial charge is found to be 

= iff + /if + /if, (4.40) 

with H\^^^ given by the restriction of the corresponding closed chain expression (see [1]), 
and 

= -{a^a^ + alal + t^t^ + rf t| + 4C/cT^Tf ). (4.41) 
Unlike the XY case, no odd-order charges exist. 
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5. Further generalizations and conclusions 

Our initial objective for tiiis work was to obtain a closed form expression for the 
conserved charges of the open XXX chain, that is, to find the finite open chain deformation 
of the Catalan tree pattern of the conserved charges for the infinite chain. This has not 
been achieved, the main reason being that the spin polynomial basis used in the description 
of the charges in the closed case is inadequate for open chains. This, however, could be 
seen only after looking at the third nontrivial conservation law (Hs). The difficulties in 
studying the open chain conservation laws also come from the absence of a boost operator; 
there is thus no recursive way of constructing the charges. 

In a sense, the problem has been solved half-way: we showed that the expressions for 
the charges split in two parts, a known bulk part and boundary corrections; we showed 
further how these corrections can be calculated systematically (and wrote down the leading 
boundary contribution) . 

Moreover, we have demonstrated that an open chain has half as many local conserved 
charges as its closed version with an identical number of sites. When one cuts open a closed 
spin chain, there appears an additional, symmetry (the chain reversal symmetry (2.22)), an 
immediate consequence of which being that all odd-order charges cease to be conserved. 
This is reminiscent of the reduction, by a factor 2, of the order of the two-dimensional 
conformal group induced by a boundary: in presence of a boundary, the holomorphic and 
antiholomorphic Virasoro modes are no longer independent! 

It should be stressed that the basic structure of the even-order conservation laws 
generated from the transfer matrix - i.e., the separation into the bulk and boundary parts 
— is a direct consequence of the construction (2.1).i The structure (2.19) characterizes 
therefore a large class of integrable chains with nontrivial boundary conditions. Moreover, 

^ As shown in sect. 2, all of the conserved charges generated by the transfer matrix of an open 
chain may be obtained in this way. Moreoer, the chain reversal symmetry guarantees that our 
construction based on solving commutativity constraints for a boundary part added to the bulk 
part of an open chain charge, cannot succeed for odd orders. Therefore, our construction provides 
no more charges then the transfer matrix formalism. However, it is not immediately clear whether 
our approach (or the transfer matrix construction) exhausts all the possible conservation laws for 
open chain. The problem of completeness of the family of conserved charges generated by the 
transfer matrix is a very difficult one and has not yet been solved for generic (closed or open) 
integrable chains; it is certainly beyond the scope of our investigation. 
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in models where the transfer matrix does not generate the complete family of conservation 
laws (in the XY model for instance), the additional conserved charges are expected to 
exhibit the same separation into bulk and boundary pieces. 

The method for the calculation of the explicit expressions of the local conserved 
charges in the open XXX model presented in sect. 3 can be applied to more general 
situations in a straightforward way. For example, through a simple unitary transforma- 
tion, the exact expressions of H4, Hq and Hg can be translated into conserved charges for 
the XXZ model with anisotropy parameter A = —1 [2]. Similarly, by reinterpreting the 
vector product and the polynomials fn,k{j) in terms of a Lie algebra commutator, these ex- 
pressions apply directly to the su{N) version of the XXX chain, as well as to the octonionic 
chain [2]. All these chains are particular examples of a general algebraic model defined 
in ref. [2], characterized by three arbitrary constants ki,K2,K3. For periodic boundary 
conditions this model has been proved to be integrable provided that ki + ks = 2k,2- The 
same condition assures also the integrability of this general model with free open boundary 
conditions. In particular, the boundary term in H4 is 



Although we have considered mainly open chains with free boundary conditions, we 
could treat more general boundary conditions in exactly the same way. Adding boundary 
terms in H2 simply induces extra terms in the higher order charges. For instance, by 
adding the boundary term kaf to the hamiltonian of the open XXX chain we find the 
following extra terms in H4: 



kalalal + kala^al - ka^a^al - kafalal - k^afa^ - k^afa^ - kaf + {Sk - k^)af. (5.2) 



Anisotropic versions of the XXX model can be studied in similar way; again, the only 
difference is that the structure of conserved charges is more complicated. For instance, the 
first nontrivial conserved charge for the XYZ chain with free open boundaries is 



with H^'^ obtained from the corresponding closed XYZ expression, (given in [1]), and 



= (2«2/«i) /2,o(l). 



(5.1) 



H^^^ = H^^^^ + h^ + h. 



R 

■4 ) 



(5.3) 



^4 = Ki^y^ + ^z^)<^l<^2 



(5.4) 
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For all these generalizations of the spin-1/2 XXX model, local conserved charges 
of odd degrees are also absent.! As already mentioned, this result is encoded in the 
transfer matrix formulation of the open chain models. Notice that when there are nontrivial 
boundary conditions, hence nontrivial matrices X±, the spin reversal transformation cr^ ^ 
aN+i-n must be accompanied by the interchange of and K_. We expect that the 
same conclusion will hold true for higher-spin chains and even nonhomogeneous chains 
with different spin representations at different sites (provide that the sequence of spaces 
^2, ■ ■ ■ 7 Vn is invariant under the interchange of sites n and N + 1 — n). 

For the simpler XY model, we could obtain the exact expression for all the charges. 
Here again, there is no open chain deformation for half of the closed chain conserved 
charges, and as for closed chains, half of the charges cannot be obtained from the power 
expansion of the transfer matrix. 

It is worth mentioning that for a particular choice of boundary terms the XXZ chain 
is endowed with quantum algebra symmetry Ug[su{2)] [13]. Clearly, the corresponding 
family of local conserved charges is then invariant under the action of this quantum algebra. 
However, this additional symmetry does not seem to lead to a significant simplification in 
the pattern of the conserved quantities. 

In ref [14], we have proposed a simple integrability test for closed chains with short- 
range intereactions. This test essentially boils down to the existence of a conserved charge 
Hs (more precisely, we conjectured that the nonexistence of Hs means nonintegrability). 
This test needs to be modified for open chains in an obvious way; here it is the existence 
of that should characterize the class of integrable models. 



6. Appendix 

For the sake of completeness, we give the basic commutators for the XXX model of 
the type [ai ■ cTi+i, fn{C)] (see eq. (4.13) in [1]) in the graphical notation of section 3. In 
the expressions below, ■ ■ ■ stand for an arbitrary number of filled (•) or open (o) dots. 

o •■■■•= 2z • (6.1) 

^ One might ask whether there are some boundary conditions such that the only charges 
present, including the hamiltonian, are of odd degree. However, this seems impossible: for ex- 
ample, for the XXX model with nonperiodic boundary conditions, the lowest order odd charges 
i^l"''^ and i/5^""= cannot be deformed in a way that preserves their mutual commutativity. 
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• o • • •• = -2z • • • • •• (6.2) 

rn 

• • •••• =4zo« •••• — 4z«o •••• (6-3) 

• ••••• •••• — 2i • ■■■o» — 2i • ••••o •••• (6-4) 

rn 

• ••••o •••• = —2i • •••(••)•••• (6-5) 



In the last equation, the parentheses indicate a modification of the usual nesting pattern, 
e.g. the sequence ••(••)••, starting at site 1, say, corresponds to the polynomial 
((((Ji X (72) X ((J3 X a^j) X (T5) ■ (Je. Further basic commutators may be obtained applying 
the chain reversal symmetry to the rules (6.1)-(6.3) and (6.5). In particular, we have: 



• . . . • o = —2i • • • • • • (6-6) 

m ■ ■ - o m = 2i • ••••• (6-7) 

• ••••• = 4i • ••••• — 4i • ■ ■ ■ • o (6.8) 

• . . . • • • • = 2z • •••(••)•••• (6.9) 
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